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A host of spatially extended systems, both in physics and in other disciplines, are well described at 
a coarse-grained scale by a Langevin equation with multiplicative-noise. Such systems may exhibit 
^ I non-equilibrium phase transitions, which can be classified into universality classes. Here we study in 

^T^. detail one of such classes that can be mapped into a Kardar-Parisi-Zhang (KPZ) interface equation 

' with a positive (negative) non-linearity in the presence of a bounding lower (upper) wall. The wall 

^ I limits the possible values taken by the height variable, introducing a lower (upper) cut-off, and induce 

a phase transition between a pinned (active) and a depinned (absorbing) phase. This transition is 
studied here using mean field and field theoretical arguments, as well as from a numerical point of 
view. Its main properties and critical features, as well as some challenging theoretical difficulties, are 
^ , reported. The differences with other multiplicative noise and bounded-KPZ universality classes are 

^>«^ ■ stressed, and the effects caused by the introduction of "attractive" walls, relevant in some physical 

' contexts, are also analyzed. 

(N ■ 

O ■ INTRODUCTION 

m 
o . 

Non-equilibrium phase transitions occurring in systems amenable to be described by Langevin equations including 
d ' a multiplicative noise (MN) term are the subject of current intense studies. This embraces a broad variety of systems 
both in physics and in other disciplines. A phenomenology much richer and complex than that appearing in equilib- 
rium systems, including counterintuitive behaviors, has been reported to appear in these, typically non-equilibrium, 
' ^ situations. See 0,0] for detailed introductions to this growing field, including many different realizations, 
g The interest in MN problems is enlarged even further, because of the existing mappings between them and other 

^ prototypical non-equilibrium problems A well known instance is the Kardar-Parisi-Zhang (KPZ) equation, de- 
IV scribing the kinetic roughening transition of generic interfaces under non-equilibrium conditions |^ , which can 
• ,—{ ', be mapped onto a MN Langevin equation, by performing the so called Cole-Hopf transformation linking the interface 
' height at each point with the activity field of the MN equation. If an interface under consideration is described 
^ , by the KPZ equation and it is physically limited by a wall, i.e. if the heights cannot be larger of smaller than a 
- - ' certain value, then this problem, bounded KPZ, can be mapped into a multiplicative noise equation by employing the 



abovementioned transformation |3j|a,l3- T'he bounded KPZ equation may experience, as parameters are varied, a 
phase transition from a depinned phase in which the interface escapes with probability one from the wall, to a pinned 
phase characterized by a finite expectation value of the stationary averaged height (measured from the wall). In the 
MN language {i.e. after employing the Cole-Hopf transformation) the pinning-depinning transition corresponds to a 
critical point separating an absorbing phase in which the order parameter goes exponentially to zero (depinned phase) 
to an active phase in which the order parameter takes a non-vanishing average value. 

Surprisingly enough, it was shown a few years ago that the introduction of "upper" or "lower" walls into a given 
KPZ equation (with a fixed non-linearity sign) lead to quite different phenomenologies. The origin of this can be 
tracked down to the fact that the KPZ equation is not invariant upon inverting the height (see U or for a more 
detailed explanation). Taking, for instance, the sign of the coefficient of the KPZ non-linearity to be positive, the 
introduction of an upper wall leads to a (well established bv now) set of critical exponents characterizing the, so 
called, multiplicative noise 1 (MNl) universality class |y, Li la |9j . On the other hand, a wall limiting negative values 
of the interface height (lower wall) leads to a different type of phase transition as shown by T. Hwa and one of us 
some years back In what follows, and following the nomenclature introduced in 01, we will use the term MN2 to 
name this class. 
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It can be easily shown that a KPZ equation with positive non-Hnearity and a lower wall is completely equivalent to 
a KPZ with a negative non-linearity coefficient and an upper wall one just have to change the sign of the height 
variable in the KPZ-like equation to verify this. 

The MN2 class is of great importance in the context of alignment of DNA and other biological sequences. It has 
been argued by Hwa and collaborators, that the phase transition appearing upon changing the, so called, scoring 
parameter in the commonly used alignment algorithm can be mapped into the MN2 critical point 0|. It is also 
related to some instances of non-equilibrium wetting For some other applications and physical instances within 
this class see and references therein. 

While the MNl class has been extensively studied, specially after its connections with non-eq uilibrium wetting 
[nl llll and with the problem of synchronization in extended systems were established 0, Q| , the MN2 class 
remains poorly studied. Furthermore, recent numerical analysis have revealed that the preliminary critical exponent 
values reported in Q might be far from their true asymptotic values. 

Aimed at clarifying these issues it is the purpose of this paper to analyze the MN2 phase transition in one dimensional 
systems using: i) mean-field and field-theoretical techniques and ii) numerical (Monte Carlo) analysis of different 
models claimed to belong to this class. 

Finally let us stress that if the wall becomes attractive, rather than simply bounding, new phenomenology might 
appear. This is particularly interesting in the context of synchronization [3.ll3l|. This possibility will also be discussed 
along the paper. 



THE MN2 CLASS 



Let us consider a KPZ equation with a positive non-linearity coefficient, A > 0, in the presence of a lower wall, 

dth{x, t)=a + b e-P'^ + DW^h + X{Wh f + ar]{x, t). (1) 

where /i is a height variable, a represents a constant drift while b e~^^ is a bounding wall. The parameter p > 
controls the wall penetrability, the limit p — s- oo corresponding to a perfectly rigid (impenetrable) wall. It has been 
shown previously that for both the MNl and the MN2 classes the magnitude of p does not influence the asymptotic 
properties at criticality 0,11 (its sign, however, is important, as it determines whether the wall is a lower or an upper 
one). The same property applies to equilibrium systems (i.e., for A = 0) 77 is a stochastic white noise with 

(r/) = and {ri{x,t)r]{x' ,t')) — 26{t — t')S{x — x'), where (•) denotes an average over the distribution of the noise. 

For a fixed value of b the interface experiences a pinning-depinning transition at some value a = ac. 

Some remarks concerning the connection of the previous equation with wetting problems follow. When h(x,t) is 
viewed as the distance separating a liquid-gas interface from a solid wall, Eq. (Q) can then be interpreted as a dynamic 
model for nonequilibrium wetting. Under this perspective a is the chemical potential difference between the liquid 
and the gas phases, (h) is the thickness of the wetting layer, and the wall is a rigid, physical substrate. At bulk phase 
coexistence, i.e. for the value of a — ac for which in the absence of the wall the interface does not move on average, 
(h) diverges at all temperatures above certain wetting temperature, Tw, while for a ^ ac the thickness of the liquid 
film can be big, but finite (pinned interface). The temperature here is controlled by the parameter 6, which vanishes 
linearly with the mean-field wetting temperature as T — Tw- Thus, on approaching coexistence for T > Tw (fo > at 
the mean-field level), (h) diverges as (h) ^ \a — a^^'^. This transition, termed complete wetting, is always continuous 
and the value of the Ph exponent depends on the nature of the forces between the particles in the fluid phases and the 
wall. In this paper only short-range, exponentially decaying interactions between all the particles and the substrate 
(as described by Eq. J^l) are considered. 

The change of variables n — exp{—h) transforms Eq. Q into the MN2 Langevin equation: 

dMx, t) = V^n - 2^ '- - (a + l)n - bnP+'^ + nr], (2) 

where for the sake of simplicity we have set A = 13 = cr = 1 and Ito calculus 16] has been used (different coefficients 
for the Laplacian and the KPZ non-linear term could be reabsorbed using n = exp(— a/i) with a proper choice of a). 
This transformation maps the depinning from the wall (h) ^ 00 to a transition into an absorbing state (n) 0. The 
physical equivalence between the cases A > with a lower- wall {p > 0), and A < with an upper wall {p < 0), reflects 
in the fact that the same equation is obtained using n ~ exp{~h) and n — exp/i, respectively. 

Observe that Eq.||2J) is identical to the MN equation describing the MNl class 0)0. Hi except for the presence of 
an extra term (Vn)^/ri. This term can also be written as (Vn) • (Vln(n)) = (Vn) • (V/i) suggesting that the interface 
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language is the natural one for this class. In fact, except for the factor 2 in front of (Vn)^/n Eq. ^ coincides with 
the Cole-Hopf transform of 

dth{x,t) ^W^h + a + be-P'' + r]{x,t) (3) 

that describes the growth of wetting layers toward their equilibrium state (observe that this is just the equilibrium, 
Edwards- Wilkinson model, in the presence of a bounding wall). Note also that the factor 2 in Eq. Q cannot be 
readsorbed by reparametrizing. 

Finally, let us underline that in the regime where 6 < the wall becomes attractive (which might be necessary to 
describe some physical situations as, for instance, synchronization problems as said in the introduction) and a new 
term, say, cexp(— 2/i) (equivalently cn^^"*"^) with c > has to be added to stabilize the equation. 

Having presented the equations defining the model, in the forthcoming sections we study the associated physics by 
using i) mean field approaches, ii) field theory, and iii) numerical, Monte Carlo simulations combined with scaling 
arguments. 



MEAN-FIELD APPROACHES 



Mean-field approaches to Eq. |^ can be implemented with several degrees of sophistication. A crude approximation 
consists of ignoring the noise and spatial variations. At this level one trivially gets that the order parameter (h) 
vanishes as a — > with an exponent — 0. When applied to Eq.Q, this approximation yields a shifted critical 
point at Oc — —1 and the usual result for the order parameter critical exponent, /3„ = 1/p. But, as experience with 
other systems with multiplicative-noise dictates, neglecting completely the noise is a too crude approximation, that 
eliminates most of the characteristic traits of MN physics. 

Now the effect of allowing a spatially varying order parameter and taking the noise into consideration are examined. 
The Laplacian is discretized as 

V^n, = ^ '^{nj - « {n) - n„ (4) 

j 

where the sum runs over the nearest neighbors of i and a large system dimensionality has been assumed. Similarly, 
the square gradient term can be written as 

A !± ^ ^ - 2{n) + n,. (5) 

The one-site stationary probability distribution is then readily obtained from the associated Fokker-Planck equation, 

Pst[n,{n)j oc ^expy ^ ^ ' — dxi, 

1 / (n) nP (n^)\ 



n 



where (n) and (n^) have to be calculated self-consistently. For a < Oc = — 5, Pst is not normalizable what means 
that the stationary state is the absorbing phase (n) = 0. For (n) ^ 0, however, the non-analyticity of exp[(n^)/n^] at 
n = again renders Pst non-normalizable. As a result, there is no well-defined active phase at this mean-field level 
for Eq. jnj. Similar problems are found if the same type of approach is applied to Eq. |^ instead of Eq. 

To avoid the presence of the square gradient term in Eq. and the complications of its, somewhat arbitrary, 
discretization, we resort to a different change of variables. After n = exp/i, 



dtn = V^n + (a + l)n + bn^^P + 
'bn^-P (n) 



2 — p 11 



nil, 

(7) 



These equations are simpler than (O and but at the cost that {n) is no longer an order parameter as, at the 
transition point (n) — s- oo rather than going to 0. Furthermore, Landau expansions only make sense when the order 
parameter vanishes at the critical point, making the whole approach inconsistent. One possible way to circumvent 
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this problem is to monitor m = 1/n, and study Pst{'m)dm = Pst{n)dn, but given the non-Gaussian nature of the 
probabiUty distribution the substitution (n) = l/(m) is hkely to be incorrect, and there seem to be no safe way to 
proceed. 

Summing up, non-trivial mean-field approaches detect some problems with the model under consideration, and are 
not able to predict a phase correct phase diagram. A sound mean-field approximation, needs therefore to be found. 
Notice that none of these problems occur in the case of a negative KPZ non-linearity, i.e. in the MNl class, where a 
standard mean- field approximation yields quahtatively correct results (see i ^ and references therein). 

FIELD THEORETICAL CONSIDERATIONS 

The Langevin equation for MNl is known to be super- renormalizable, i.e. Feynman diagrams can be computed 
to all orders and resummated. This does not imply that critical exponents can be computed in all the cases, as the 
renormalization group flow-equation has runaway trajectories supposed to converge to a .strong coupling fixed point. 
But at least, the correct phase diagram, including strong and weak coupling fixed points can be obtained. 

For the MN2 the situation is far more complicated, as can be a priori anticipated given the failure of mean-field 
approaches. The extra term, (Vn)^/n, being singular in n precludes the use of perturbative expansions around n — 0. 
Given the lack of a non-pcrturbative approach to the KPZ and MN strong coupling fixed points, there is not much we 
can add to this section, except that there is a promising attempt to tackle this and related KPZ-like problems. There 
is a formalism, developed by Fogedby, aimed at developing a strong coupling theory for KPZ based on a semiclassical 
or WKB approximation applied upon the Martin-Siggia-Rose generating functional 20]. Its main advantage is that 
it does not involve expansions around classical noiseless solutions, but around classical (extremal) noisy solutions. It 
would be very interesting to extend these ideas to KPZ problems in the presence walls, namely, to the multiplicative 
noise universality classes MNl and MN2. 

NUMERICAL RESULTS 

Owing to the failure of standard (non-trivial) mean-field approaches and lacking so far of an alternative analytical 
route, numerical methods are required to glean insight into the system properties. We have carried out simulations 
of two surface growth models. Both of them, in the absence of walls, are known to belong to the KPZ universality 
class. An extra rule is then added to generate a bounding wall, as described above. 

Model 1 

Our first model was introduced in and is defined as follows: the surface position at time t above a site x on a 
one-dimensional lattice of size L is given by a continuous height variable tit{x). A new height configuration is then 
generated in a three-step process. 

1. Each lattice site ht{x) is updated according to h^{x) = ht{x) + a + rjtix). 7]t{x) is a random number uniformly 
distributed in [0,1] and a is a constant drift term analogous to that of l[T|l. 

2. The configuration is changed to ht+i{x) — min[/i((a: ± 1) -I- 7, /i((a;)], where 7 is a constant whose precise value 
is not essential for the final results. We have set without loss of generality 7 = 0.1 as in previous numerical 
analyses. 

3. A hard wall at = is introduced by way of the additional rule }it{x) ~ min[ft,i(a;), 0] |^. 

Finally, periodic boundary conditions are imposed and tit{x) is initially set to 0. 

The continuum counterpart of this model is known to be a KPZ equation |2ll ] with A < in the presence of an 
upper wall '22] which, as remarked before, is equivalent to the case A > and a lower wall, and corresponds to the 
MN2 class. Numerical results for this model were first presented in Q and seemed to be consistent with mean-field 
{i.e. single-site) like behavior. However, a similar problem recently studied in the context of synchronization has 
revealed inconsistencies probably due to insufficient statistics 0, ll^ . In this subsection improved simulation results 
are provided upon revisiting the analysis reported in 8] for larger system sizes and longer sampling times. The case 
of an attractive wall, not included in is also considered. 
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FIG. 1: Model 1: steady-state values of the oder parameter, (n) = (exp/i), as function of the system size, L, for drifts (top 
to bottom) 1.57450, 1.57440, 1.57435, 1.57433, 1.57430, 1.57425, 1.57420. (•) denotes both spatial and temporal averages, as 
well as averages over independent runs. The straight line corresponds to the critical point = 1.57433(2) and from its slope 
f3„/iy± ^ 0.33(2). 

First, we take up the case of a simple (non-attractive) wall, corresponding to 6 > in Eqs. ^ and Q. Figure (1) 
shows how the steady order parameter (n) changes with the system size L. Within the active phase, it saturates to a 
constant value, while in the absorbing one it bends down and decays exponentially. Our best estimate for the critical 
point is Qc = 1.57433(2) and from the slope of the curve we get l3n/i^± ~ 0.33(2) This value of ac corresponds 
to the point where a free interface (far from the wall) has zero average velocity. The time evolution of (n) at the 
critical point for different system sizes (Figure (2)) behaves like (n(t)) ~ with 9n = 0.215(15) or, equivalently, 

{h{t)) ~ t^*^, with 6h ~ 0.355(15). As for the exponents /3, which govern the saturation of the order parameter within 
the active phase, they have been computed using the largest available system sizes {L = 1600 and 3200). The best 
fit to {X) ~ |a — ad^^^ yields /3„ « 0.32(3) (for (n)) the and Ph = 0.52(2) (for (h)), respectively. The error margin 
is typically larger here than for other exponents due to the sensitivity to the uncertainty in the determination of the 
critical point. 

It was proved in [^Q that these exponents must satisfy the scaling relation z — (3/{v^6), where z is the dynamic 
exponent of the KPZ equation. The exact value for zvad— 1 is 3/2, thereby z = 0.33/0.215 = 1.5(1) in agreement with 
the prediction. In addition, also from [^0, i^j. = 1 , and from our direct measurements we get [Pn/ ^i.) / Pn = 0.33/0.32, 
implying v± « 1. In terms of h and assuming v± = 1^ z — 0.52/. 355 = 1.5(1) which, again, is compatible with 3/2 
within error bars. 

The two alternative, but equivalent, mathematical descriptions of the MN2 class in terms of h and n = exp(— /i) 
can be related noting that the latter is essentially the density of sites at zero hight, n{x,t) = (5/i(2;.t),o 24]. We have 
verified that n and 5i^i^x t).o exhibit the same asymptotic scaling behavior. 

Attractive walls can also be simulated within this model by simply substituting a by a — bSh.o, where 6 < and the 
sign convention is chosen to keep the analogy with Eq. This means that wherever the interface in attached to 
the wall it experiences an additional ("sticky") force pushing it against the wall. Extensive Monte Carlo simulations 
for b = —0.3 show that the previous results, in what respect universal features, carry over without change, the 
only difference being that the approach to asymptotics is slower. Upon increasing the attractiveness of the wall the 
transients become longer. The estimate for the critical point is the same one as before (this is due to the fact that the 
free interface is not affected by variations of the attractiveness parameter). Again, our best estimates for the critical 
exponents are (in/^A. — 0.32 and 6'„ — 0.215 (see Figure (3)). For b = —0.4 we still observe a second-order phase 
transition with a crossover to the mentioned exponents. For b = —0.5 transition becomes first-order but it still occurs 
at a = ac. Within the active phase all the sites are closed to the wall and the order parameter is 1, but it suddenly 
changes to upon decreasing a and hysteresis is observed for slightly subcritical values a (the interface is pinned for 
up to long times). Therefore, a tricritical point must exist between 0.4 and 0.5. We have identified it at 6 = —0.42(1). 
The tricritical behavior has not been analyzed. 

Let us stress that, contrarily to what happens for the MNl class, where the presence of an attractive wall induces a 
new and rich phenomenology (including a broad region of phase coexistence and directed-percolation type of transitions 
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FIG. 2: Model 1: time evolution of the order parameter, {n{t)), at criticality, a = Uc, for system sizes 100, 400, 800, and 1600. 
The saturation values are those of Figure 1. The best fit gives 9n = 0.215(15). The straight line is a guide for the eye and has 
a slope -0.215. 




FIG. 3: Model 1: steady-state values of the oder parameter for an attractive wall, b = —0.3, as function of the system size 
L for drifts (top to bottom) 1.57450, 1.57447, 1.57445, 1.57440, 1.57435, 1.57433, 1.57430, 1.57425, 1.57420. The best fit to a 
straight line again corresponds to Oc = 1.57433. 

), the addition of "attractiveness" has a very mild effect here. Basically, it just shifts the position of the critical 
point and induces a first-order phase transition for very strong attractions. 

Model 2 

In order to verify the robustness and eventual universality of the previous results, we have performed a second study 
of a different model. It is a restricted solid-on-solid (RSOS) model, a variant of the single-step model introduced in 
^ , in the presence of a wall. A similar model has been recently studied in the context of synchronization transitions 
25| . to study MNl type of transitions. Initially the wall is located at = and a grooved interface is placed beneath 
it, i.e. the interface has negative height at all the positions. 

The dynamics proceeds as follows: At each time step, a site is randomly picked from a one-dimensional lattice of 
length L and its height decreased two units, h{i) — > h{i) — 2, provided that h{i) > h{i + 1) and h{i) > h{i — 1), 
i.e. provided that it is a local maximum. Should the the RSOS constraint be violated, the trial is discarded and 
repeated. Every 2(L — 1)/[1 — 25v{l — steps the wall retreats one unit and, simultaneously, the interface is moved 
downwards by two units wherever it lies above the wall |25l| . The difference between the wall and interface velocities. 
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Model 1 0.32(2) 0.97(5) 0.34(2) 1.55(5) 0.215(15) not measured 
Model 2 0.325(5) ^ 1 0.33(2) ki 1.5 0.215(5) rj = 0.8 

TABLE L Critical exponents for the MN2 class in d = 1. 



A B 
0| 1 1 1 Ip™ r—- 




FIG. 4: Model 2; (A) order parameter behavior, {exp{hw — h)), in the vicinity of the critical point Svc — 0. From the slope of 
the line we get /3„ = 0.325(5). (B) Decay of the order parameter, {exp{hw ~ h)), at the critical point yields 6n = 0.215(5) (cf. 
with Figure (2). The straight line is a guide for the eye and has a slope —0.215. 

Sv, acts as the control parameter: if 6v is negative, the interface eventually depins from the wall, while for 6v > it 
remains pinned). It can be easily shown, using random walk arguments that for the chosen wall velocity the system 
seats at its critical point. By varying it, we have a control parameter. The possibility of computing analytically the 
critical point largely simplifies the numerical analysis, and makes of this an efhcient discrete model. 

The quantities monitored are (exp(/ii„ — h)) and /i^ — h. We have measured the exponents /?„ and 0„ for a system 
of L = 2^° sites. Our results lead to /?„ = 0.325(5) and 6'„ = 0.215(5), in excellent agreement with those of model 1 
(Figure (4)). Once more, assuming i^^ ~ 1, we get z k, 1.51 in good agreement with the scaling laws. In addition, we 
have also measured the spreading exponent, 77, that characterizes the number of pinned sites. It is computed averaging 
over all the runs and starting with an initial condition with a single point attached to the wall. Our best estimate 
is 77 = 0.80(2) (Figure (5)). Measuring the surviving probability, and therefore the exponents 8 and CJ , is a delicate 
technical point because it is hard to decide when the activity of a run has ceased. We have not tackled this problem 
here. Lastly, we obtain — 0.34(5), which is again in good agreement with the value reported for our Model 1. 

We have also considered different variations of the model in which the interface can penetrate the wall at some 
points, i.e. the wall is not perfectly rigid. None of the universal features seem to be affected by this change. 

In conclusion, all these results support strongly the existence of robust universality in the MN2 class. 

As a matter of consistency, we have modified the algorithm of the model to simulate a lower wall, and therefore 
a case expected to be in the MNl class. We obtain the set of exponents /?„ — 1.69, = 1.19 and r] = —0.4, all of 
them in agreement we Dreviously reported results and showing that the upper and lower problem belong to different 
universality classes 0, B 1^] • 

Numeric integration of stochastic differential equations 

As a further test for universality, we have numerically integrated Eq. using a Milshtein's algorithm ^26] . A system 
size of L = 2000 was considered and the time step and mesh size were set to At = 0.001 and Ax = 1, respectively. For 
simulation times up to i = 10^ (10^ trials per site) our results lie far from the asymptotic regime. We do not discard 
numerical instabilities in the integration scheme, as it is known that results from numerical integrations may not agree 
with the predictions from the continuum KPZ QJ^. The Cole-Hopf transform is a standard way to account numerically 
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FIG. 5: Model 2: number of pinned sites as a function of time for an initial condition with a single point attached to the wall. 
A fit for late times yields rj = 0.80(2). The straight line is a guide for the eye and has slope 0.8. 

for the integration of bounded KPZ equations and, indeed, Langevin equations with MN are by far more stable than 
their interface-language KPZ-like counterparts 0, Nevertheless, we have also found numerical problems when 
integrating either when the extra term is written in logarithmic form or averaging for smoothing the gradient. 
All numerical attempts are unstable nearby the absorbing state, owing to the presence of the extra singular term. We 
leave, therefore, the numerical integration of a continuous Langevin equation, representative of the MN2 class as an 
open, challenging problem. 

DISCUSSION 

We have characterized the MN2 universality class, or analogously its bounded KPZ counterpart, which, as com- 
mented above, accommodates different physical phenomena. We have studied it from, somehow deceptive mean-field 
and field theoretical approaches, as well as by numerical studies. None of the analytical methods provides a satis- 
factory description of the phase transition present in this class. On the other hand, Monte Carlo simulations of two 
different discrete interface models, argued to belong to this universality class, give a firm evidence for the existence of 
a robust universality class. Contrarily to what previous simulations seemed to indicate our results are not a simple 
extension of the ones obtained for one-site, implying that spatial correlations play an important role. Table 1 gathers 
the values of the critical exponents in terms of n for the two discrete models considered in this paper. From the Monte 
Carlo estimates, it cannot be discarded that they adopt the rational values /?„ = 1/3 and 6^ ~ 2/9, which combined 
with the exact values derived in i'± — 1 and z = 3/2, would lead to Bn/ vi. = 1/3. Note that our results do not 

compare well with those of the nonequilibrium wetting model reported in . We believe this is probably due to the 
extremely long transients known to be present in that model. We have verified that the transition point is located 
at the same value of the control parameter for any value of the "attractiveness" parameter (6 in Model 1), either 
representing an attractive or a non attractive wall. For strong enough attractive walls, i.e. b sufRciently negative, the 
transition becomes first order as in 11], while if the wall is weakly attractive then it remains in the MN2 class. 

The problem of reaching a satisfactory analytical understanding, and even that of obtaining sound results from 
numerical integrations of the continuous Langevin-equation (in either the interface or the density language) represen- 
tative of this class remains an open challenge. 

Summing up, even though strong evidence is provided confirming the existence of a universality class (i.e. the 
corresponding critical exponents are computed with good precision in one dimension and they arc universal in two 
different discrete models), its theoretical description in terms of Langevin equations, contrarily to what happens for 
the closely related MNl class, is far from satisfactory. In particular, the Langevin equation does not seem to admit 
sound mean-field solutions, nor is amenable to be treated by means of standard perturbative field theoretical tools, nor 
it admits a stable numerical integration. Identifying the physical causes at the root of these difficulties is a challenge 
for future research. 
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